Exact diagonalization study of optical conductivity in two-dimensional Hubbard model 
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The optical conductivity cr(u>) in the two-dimensional Hubbard model is examined by applying 
the exact diagonalization technique to small square clusters with periodic boundary conditions up 
to v20 x \/20 sites. Spectral-weight distributions at half filling and their doping dependence in the 
20-site cluster are found to be similar to those in a \/l8 x \/l8 cluster, but different from 4x4 results. 
The results for the 20-site cluster enable us to perform a systematic study of the doping dependence 
of the spectral-weight transfer from the region of the Mott-gap excitation to lower-energy regions. 
We discuss the dependence of the Drude weight and the effective carrier number on the electron 
density at a large on-site Coulomb interaction. 

PACS numbers: 71.10.Fd, 78.20.Bh, 78.20.Ci 



I. INTRODUCTION 

Charge carriers introduced into the Mott insulators in- 
duce a dramatic change of the electronic states. The 
high-temperature cuprate superconductor is a good ex- 
ample, where the spectral weight of the Mott-gap exci- 
tation around ~ 2 eV in the optical conductivity, u{oS), 
decreases rapidly upon carrier doping and new weight 
emerges in the mid-infrared region accompanied with the 
increase of a zero-frequency component that corresponds 
to the dc conductivityi2i2ii Such a reconstruction of the 
electronic states would not occur in a conventional semi- 
conductor. 

The electronic states of the cuprates are known to be 
well described by a single-band Hubbard model with a 
large on-site Coulomb interaction in two dimensions. Nu- 
merous theoretical works have been done to understand 
the physics of the two-dimensional (2D) Hubbard model. 
Among them, the numerical diagonalization technique 
based on the Lanczos algorithms has been intensively 
used for the investigation of the rapid spectral weight 
transfer seen in <j(lo) and other quantities^Si&Si The 
doping dependence of <j{u>) for a 4 x 4 cluster under 
periodic boundary conditions has clearly shown the re- 
construction of the weight^ However, the cluster signifi- 
cantly suffers from the finite-size effect. For example, the 
Drude weight D at half-filling, where the electron density 
n = 1 (16 electrons in the cluster), becomes negative, and 
the ground state at n = 0.875 (14 electrons) is unphysi- 
cally threefold degenerated. Although D seems to grow 
proportional to 1 — n near half filling^ this would not 
be conclusive because of these finite-size effects. It has 
been argued^ that D in a 4 x 4 cluster with mixed bound- 
ary conditions rather grows proportional to (1 — n) 2 . In 
these situations, it is necessary to perform a systematic 
study of the optical conductivity by increasing the size 
of clusters. 

In this paper, we perform the exact diagonalization 
study of <r(w) in the 2D Hubbard model up to a 20-site 
cluster. The diagonalization of the 20-site cluster has not 
been done in the past as far as we know. We compare 
a{uj) for 4 x 4, \/l8 x y/18, and y/20 x \/20 clusters with 



periodic boundary conditions. We find that at half-filling 
the 18- and 20-site results are similar in spectral shape 
but different from the 16-site result about the spectral 
weight distribution. The doping dependence of a(u>) also 
shows similar behavior between the 18- and 20-site lat- 
tices. For the 20-site cluster, we perform a systematic 
study of the doping dependence of the spectral-weight 
transfer, and find that the transferred weight from the 
Mott-gap excitation increases with doping, showing a be- 
havior qualitatively consistent with that observed exper- 
imentally- We also find that D grows roughly propor- 
tional to 1 — n near half-filling. These exact results can 
be useful in understanding the physics of the 2D Hubbard 
model. 

This paper is organized as follows. In Sec.|H] we intro- 
duce the 2D Hubbard model and the method for the exact 
diagonalization. The doping dependence of the ground- 
state energy is also presented in the section. In Sec. IIIII 
the calculated results of a{uS) are shown. The doping de- 
pendence of the spectral- weight transfer, as well as of the 
Drude weight, is performed for the 20-site cluster. The 
summary is given in Sec. IIVI 

II. MODEL AND METHOD 

The Hamiltonian of the Hubbard model in 2D reads as 

i,S,a i 

W 

where t is the hopping energy, U is the on-site Coulomb 
interation, and S = x and y, x and y are the vec- 
tors connecting the nearest-neighbor sites in the x and 
y directions, respectively. The operator c\ a (c- 1:a ) cre- 
ates (annihilates) an electron with spin a at site i, and 

We apply the exact diagonalization method based on 
the Lanczos technique to square clusters with periodic 
boundary conditions up to 20 sites. The size of the 
Hilbert space to be diagonalized can be reduced by em- 
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FIG. 1: Ground-state energy Eq (a) and energy of the hop- 
ping part (kinetic energy) K (b) per site for the %/20 x \/20 
cluster (AT = 20) with an even number of electrons. U/t — 5 
(circles), 10 (upper triangles), and 20 (squares). The energies 
for U/t — (solid lines) are obtained from the infinite-size 
system. 



ploying the translational and point-group symmetries of 
the clusters. In a 4 x 4 cluster, the maximum size of 
the Hamiltonian matrix is 5,178,144 for the space with 
n — 1 (half-filling), null z-component of the total spin 
(S z = 0), the total momentum K = (±tt/2, 0), (0, ±tt/2), 
(it, ±7t/2), and (±7r/2, tt) with Ai irreducible representa- 
tion. Since the memory to be required for a single vector 
is 79 MB (complex with double precision), a standard 
workstation (Pentium 4, 2.67 GHz, 1 GB memory) is 
enough for the Lanczos diagonalization: the ground-state 
energy is obtained in 4 h by using our code. The size of 
the matrix rapidly increases with the size of cluster. In 
a V20 x V20 cluster, the maximum size is 1,706,748,500 
for the space with n = 1, S z = 0, K = (±47r/5, ±27r/5) 
and (±27t/5,=f47t/5) (25.43 GB for a single vector). For 
such a huge matrix, we made use of the supercomputer 
(Hitachi SRI 1000 in the Institute for Molecular Science, 
Okazaki) and a parallelization based on the message pass- 
ing interface. In order to obtain the ground-state energy, 
it takes 37 h when we use 16 nodes in the supercomputer. 

The ground-state energy E of the y20 x y20 cluster 
with an even number of electrons is plotted in Fig. ^a) 
as a function of n for U/t — 5, 10, and 20. In the weak- 
coupling regime with a small value of U, the ground- 
state energy can be obtained by ordinary perturbation 
theoryJ^. For U/t >5 we find no indication of phase 
separation from Fig. ^a). We also find that the energy 
curve does not behave smoothly at n = 0.5. This devia- 
tion probably comes from a closed-shell structure formed 
in the momentum space when n = 0.5 (ten electrons in 
the cluster). In other words, the kinetic-energy term may 
be responsible for the deviation. In fact, the same devi- 
ation is seen in the kinetic energy, as shown in Fig.^b). 

As a database, we have listed in Ref. 11 the ground- 
state energies of the \/l8 x yl8 and \/20 x \/20 clusters 
with U/t — 5, 10, and 20 for all n, together with the 
lowest eigenenergy in each subspace of the clusters. We 
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FIG. 2: Regular part of the optical conductivity a Icg (ui) at 
half-filling (n = 1). (a)-(c) U/t = 5. (d)-(f) U/t = 10. (a) 
and (d) 4x4. (b) and (e) 718 x v^IS. (c) and (f) \/20 x ^20- 



note that for U/t = 5, 10, and 20 the ground states of the 
V20 x \/20 cluster with an even number of electrons are 
always located in the subspace of zero total momentum, 
but not for v/l8 x v/l8. 



III. OPTICAL CONDUCTIVITY 

The real part of the optical conductivity under the 
electric field applied along the x direction is given by 
<j(u>) = 2ttDS(uj) + a Teg (u)), where D is the charge stiff- 
ness, which is sometimes called the Drude weight. The 
regular part <7 reg (a;) is 

a Ies (co) = -^-J2 \(^m\j x \^o)\ 2 S(io~E m + E ) , (2) 



where N is the number of sites, (^m) is the ground 
(final) state with energy Eq (E m ), and the electric charge 
e, the lattice spacing a, and h are taken to be unity. 
The x component of the current operator is given by 



rule, D is given by 



r Ci iCr ). From the /-sum 



D 



1 

~2N 



a leg (uj)duj , (3) 



where the stress-tensor operator t xx is given by t xx — 

"^^(Ci+x^^^+'L/i,")- We note tnat tne double 
of the first term on the right-hand side of Eq. © is plot- 
ted in Fig. njb) f° r the V20 x V20 cluster. A standard 
Lanczos technique with a Lorentzian broadening of 0.2£ 
is used in calculating Eq. @. 

Figure|2shows the dependence of cr r eg(w) on the size of 
clusters at half-filling with U/t = 5 and 10. Irrespective 
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of the size, there is a sharp peak at the Mott-gap edge. 
This is an cxcitonic peak caused by the spin-polaron for- 
mation in the photoexcited state i 12 ' 1 ? We note that such 
an excitonic peak is not seen in cr re g(w) of the ID Hub- 
bard model at half-filling. 14 We find in the figure that the 
weight of the peak is dependent on the system size: The 
weights for VT8 x \/l8 and \/20 x \/20 are smaller than 
that for 4x4. Size dependence is also seen at a higher- 
energy region, where broad spectral weights spread: The 
4x4 result exhibits a pronounced peak structure, for ex- 
ample, at oj/t = 9 for U/t = 10, which is absent in the 
\/T8 x a/18 and V20 x V20 cases. From the figure, we 
can conclude that the a/18 x ^/l8 and V20 x V20 results 
are qualitatively the same as that of 4 x 4, but quan- 
titatively different. Such a quantitative difference may 
partly come from a special shape of the 4x4 periodic 
cluster that possesses the symmetry of a hypercube^ 

The doping dependence of a reg (u) near half-filling for 
VT8 x ^/T8 and V20 x v 7 ^) are shown in Fig. 03 For all 
values of U/t shown in the figure, the weight of the Mott- 
gap excitation decreases by decreasing the density n from 
n = 1, being consistent with the previous 4x4 resultsi££ 
The most pronounced feature is the disappearance of the 
excitonic peak seen at half-filling. This is understood as a 
result of the suppression of antiferromagnetic order, since 
the excitonic-peak formation is closely related to the spin 
correlation in the spin backgrounds^ The suppression 
of the peak with doping has recently been observed in 
Na-doped Ca 2 Cu0 2 Cl 2 ^ 

In contrast to the suppression of the Mott-gap exci- 
tation, the spectral weight emerges inside the gap with 
carrier doping. The weight gradually increases with de- 
creasing n in both the V18x \/l8 and \/20x \/20 clusters. 
This is clearly seen in the cases of U/t = 10 and 20 be- 
cause of an obvious separation between the in-gap and 
Mott-gap excitations. For U/t = 5, the separation is not 
clear, however at around oj/t = 2 the weight increases 
with decreasing n. We note that in the ID Hubbard 
model, no spectral weight emerges inside the gap upon 
carrier doping^ 

In a 4 x 4 cluster with periodic boundary conditions, 5 
the weight inside the gap was reported to increase with 
reducing n, but the result was obtained by subtracting a 
lowest-energy peak with large weight that was regarded 
as a part of the Drude weight. On the other hand, the re- 
sults for a 4 x 4 cluster with mixed boundary conditions^ 
showed a decrease of the in-gap weight with decreasing 
n. This contradictory behavior in the two 4x4 clus- 
ters seems to come from the finite-size effect that is en- 
hanced in the case of 14 electrons in the clusters. In our 
Vl8 x \/l8 and -s/20 x V20 clusters, however, the doping 
dependence of the in-gap weight near half-filling shows 
the same tendency and is consistent with the experimen- 
tal dataiSi&i We note that a dramatic increase of the 
low-energy peak at n = 12/18 ~ 0.667 in Fig. ^c) is 
probably due to nonzero total momentum of the ground 
state [K = (±27r/3,±27r/3) (Ref. 11)]. 

In order to see the spectral-weight transfer more 
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FIG. 3: (Color online) Dependence of the regular part of the 
optical conductivity a rcg (cu) on the electron density n near 
half-filling, (a)-(c) Vl8 X Vl8. (d)-(f) y/20 x V%0. (a) and 
(d) U/t = 5. (b) and (e) U/t = 10. (c) and (e) U/t = 20. The 
black (dash-dotted), green (gray dash-dotted), blue (solid), 
and red (dark-gray solid) lines in the 20-site (18-site) cluster 
correspond to the results of n = 1 (1), 0.9 (0.889), 0.8 (0.778), 
and 0.7 (0.667), respectively. 



clearly, we plot in Fig. 0] integrated weights above and 
below the Mott gap for the \/20 x \/20 cluster. The in- 
tegrated weight above the Mott gap, Wm, is obtained 
by subtracting the integrated weight below the Mott-gap 
excitation, Wl, from the total weight of a reg (uj). The 
Wm decreases rapidly by decreasing n from half-filling, 
and almost disappears below n = 0.3. This behavior 
is consistent with that of an analytical expression that 
has been derived from the large U limit of the Hubbard 
model^ where the nearest-neighbor spin-spin correlation 
as well as three-site hopping processes of electrons pre- 
dominantly determines Wm- The Wl, on the other hand, 
increases with decreasing n from n = 1, and shows two 
maxima at around n = 0.7 and 0.5. The latter max- 
imum might be related to the closed-shell structure in 
the case of ten electrons, as mentioned above. Regarding 
the n = 0.5 maximum as a consequence of the finite-size 
effect, we find that the doping dependence of Wl is qual- 
itatively similar to that in the infinite-dimensional Hub- 
bard models where Wl show a broad peak at around 
n = 0.75. 

The Wl is enhanced with increasing U, i.e., with de- 
creasing the AF exchange constant J = 4t 2 /U, as seen in 
Fig. This enhancement is accompanied by a reduction 
of D, as will be shown below. The decrease of J may in- 
duce the suppression of coherent motion of spin-polaronic 
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FIG. 4: Dependence of integrated spectral weights in a Tes (u>) 
on the electron density n for the a/20 x \/20 cluster. The 
weights of the Mott-gap excitation, Wm, are shown by the 
open upper triangles and squares for U/t — 10 and 20, re- 
spectively. The solid upper triangles, squares, and diamonds 
represent the low-energy weight below the Mott gap, Wl, for 
U/t = 10, 20, and 100, respectively. 
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FIG. 5: Dependence of the Drude weight D (a) and effective 
carrier number N e g below the Mott gap (b) on the electron 
density n for several values of U/t in the \/20 x \/20 cluster. 
U/t = 5 (circles), 8 (lower triangles), 10 (upper triangles), and 
20 (squares). The solid line is for U/t = 0. The open circles in 
(b) correspond to the experimental data for L^-ajSr^CuO^ 
(Ref. 1). 



particles but enhance their incoherent motioniii This re- 
sults in the enhancement of Wl together with the reduc- 
tion of D. 

Figure [Sfa) shows the n dependence of D for the 20- 
site cluster. For large U/t, D at half-filling (n = 1) is 
close to zero, in contrast to the case of the periodic 4x4 
cluster, where D is negative at n = 1. Therefore, the 
finite-size effect on D is remarkably reduced. However, 
in U/t = 5, D at n = 1 is finite with large value, being 
different from a plausible view of an insulating ground 
state at half-filling. 18-19 This implies that the finite-size 
effect remains large for a small value of U. For very 
large U, D exhibits a dome-like shape as a function of 
n. The maximum position of D shifts toward n — 1 with 
decreasing U, as expected. 

The doping dependence of D near half-filling would be 



related to the nature of the metal- insulator transition. 20 
It has been argued 6 that, at U/t = 16, D in the 4x4 clus- 
ter with mixed boundary conditions grows proportional 
to (1 — n) 2 near half-filling by examining the values of 
D at n = 1, 0.875, and 0.75. In contrast to this, our 
data for large U in Fig. Eta) grow roughly proportional 
to 1 — rt, with a slight convex behavior from n = 1 to 
0.8. , i.e., D oc (1 - n) a with a sC 1. The results of the 
\/l8 x V4~8 cluster show the same behavior (not shown) . 
In the strong coupling limit of the 20-site Hubbard model 
(the t-J model plus three-site terms), the tendency is un- 
changed under mixed boundary conditions* 21 A similar 
dependence has also been reported in a dynamical mean- 
field theory to incommensurate magnetic ordered states 
in the 2D Hubbard model. 22 However, since our clusters 
are still small, we need to work on larger systems to ob- 
tain more conclusive results about the n dependence of 
D near the metal-insulator transition. 

In Fig. GJb), we plot the effective carrier density N c g 
for the a/20 x \/20 cluster, defined as£ 

N cs = ^ + a leg (uj)duij . (4) 

where uj c is an energy just below the Mott-gap excita- 
tion. The dependence of N c g on n shows a behavior 
qualitatively similar to that observed experimentally for 
L^-zSr^CuC^i N e ff has a tendency to be saturated at 
large doping. The experimental data seem to be fitted 
well to the U/t = 8 results, but for a more complete com- 
parison we need to include, for example, the long-range 
hoppings that are necessary for the study of the material 
dependence^ 



IV. SUMMARY 

In summary, we have examined ct(uj) in the 2D Hub- 
bard model by applying the exact diagonalization tech- 
nique to small square clusters up to 20 sites under pe- 
riodic boundary conditions. The 20 sites in 2D are the 
largest size treated ever by the exact diagonalization as 
far as we know. We have found that the spectral shapes 
of cr(cj) at half-filling and their doping dependence in the 
20-site cluster are similar to those in the 18-site cluster, 
but different from the 16-site results. By using the 20- 
site results, we have analyzed the doping dependence of 
the spectral- weight transfer as well as of D, and clarified 
systematic changes of <t(uj) in the 2D Hubbard model, 
which provide important information on the physics of 
the model. 
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